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ABSTRACT 

We point out that already existing literature on relativistic coUisionless 
MHD shocks show that the parameter a = upstream proper magnetic 
energy density/upstream rest mass energy density, plays an important role 
in determining the structure and accelerating properties of such shocks. By 
adopting the value of cr fa 0.002 which corresponds to the relativistic shock 
associated with the Crab nebula, and by using appropriate relativistic shock 
jump conditions, we obtain here a generous upper-limit on the value of (proper) 
magnetic field, ~ 1-5 x 10~^77n}^^G, for gamma ray burst (GRB) blast wave. 
Here, rj = E/Mc^, where E is the energy and M is the mass of the baryons 
entrained in the original fireball (FB), and rii is the proper number density of 
the ambient medium. Further, wc point out that, in realistic cases, the actual 
value Bsh could be as low as ~ 5 x lQ~^7]ni ^G. realistic cases. 

Subject headings: gamma rays: bursts- hydrodynamics-relativity - shock waves 
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1. Introduction 

Understanding the phenomenon of GRBs is one of the important problems of recent 
astrophysics. Fortunately, following the discovery of a cosmological redshift in the May 
08, 1997 event, it is certain now that some or all of them could be of cosmological origin 
Metzger et al. (1997). Whether cosmological or galactic, GRB phenomenon is broadly 
understood in terms of a standard model developed by Cavallo & Rees (1978), Goodman 
(1986), Paczynski (1986), Eichler et al. (1989), Shemi & Piran (1990). Nonetheless, as far 
as the origin of the complex nonthermal observed GRB spectra are concerned, an important 
development took place with the work of Rees & Meszaros (1992) and Meszaros & Rees 
(1993) suggesting that the cosmic fireballs (FB) with an optimal amount of baryonic 
pollution, T] ~ 10^ — 10^, could explain such spectra, where r] = E/Mc^, where E is the 
energy and M is the mass of the baryons entrained in the original fireball (FB). Meszaros 
& Rees suggested that as the baryon polluted FB deposits half of its original momentum 
onto the ambient medium, presumably, the bare interstellar medium (ISM), at r = r^, the 
so-called deceleration radius, the blast wave becomes very strong and radiates part of its 
energy. For further appreciation of this paper it would be appropriate to crudely visualize 
the geometry associated with the blast wave in terms of a 1-D simple diagram (Mitra 1998, 
henceforth M98). Here region (1) is the ambient ISM, the lab frame, 5*1 is the forward shock 
moving ahead of the contact discontinuity S, the location of the original FB boundary. 
Region (2) represents the (forward) shocked fluid and it is this region which is the site for 
the particle acceleration and gamma ray production in this standard model. The region (4) 
is the unperurbed FB and (3) is the part of the FB compressed by the reverse shock front 
5*2. It was shown in M98 that, in the context of this standard model, the reverse shock 
plays an insignificant role in the overall energy balance and may be neglected for dynamical 
purposes. 
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The gamma rays are likely to be produced either by a synchrotron process or a 
self-synchrotron-Compton process occuring near the region (2) and the most crucial factor 
for the success of such processes is the value of the comoving magnetic field B'2 = Bgh- 
Here prime denotes respective comoving quantities, i.e. respective proper values, whereas 
'*' denotes quantities measured in the rest frame of the forward shock 5"!. The question of 
probable generation of a magnetic field in a relativistic (or nonrelativistic) shock is a poorly 
understood topic, and, practically, most of the authors have therefore been compelled to 
use an equiparition argument to estimate the same (Meszaros & Rees 1993, Cheng & Wei 
1996, Vietri 1995, Waxman 1995): 



1^ - V'n.mc' (1) 



But by recalling the basic definition of 



3E \ 



7 X IQ^'Eli^fW cm (2) 



where rj^ — 10~^rj and £'51 = £'/10^^erg s"^, it can be easily verified that the energy density 
shown on the R.H.S. of eq.(l) directly corresponds to the region (4), i.e, the unperturbed 
FB, having a proper density (M98) 



< = ^ 5 X lQ'E,,rii{^^/v) cm-^ (3) 



Anr^c^m 



i.e., actually. 



^ 7|nimc2 (4) 



Here 7f is the bulk Lorentz factor (LF) of Si in the lab frame (1), m is the mass of a 
proton, and rii is the particle number density of the ambient medium in units of lcm~^. 
Note that for the region (1) the comoving frame coincides with the lab frame and at r = ra, 
we will have 7f ~ 77/2 in the Meszaros & Rees scenario. On the other hand, we need to 
apply the equipartition argument in the downstream of the shock, i.e., in the region, which 
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is expected to be turbulent, and which is, in any case, the site for the particle acceleration. 
And, it will be seen in the next section that, cquipartition argument yields approximately 
the same value of B2 in the shocked fluid, and this is important, because, energy density, in 
itself, is not a Lorentz invariant quantity. 

Nevertheless, the question we want to pose here is how justified is this assumption 
of cquipartition in the context of GRBs and whether by adopting this brute assumption 
we are running into conflict with some well established feature of relativistic coUisionless 
shocks. As was stressed in Mitra (1996), cquipartition, as a general physical concept may 
be found to be valid in steady-state situations like the ISM where the plasma interacts 
with the particles and currents over astronomically significant time scales. As to dynamic 
situations, there are hints that many young supernova remnants are endowded with 
freshly generated magnetic fields which are considerably higher than the bare ISM values 
~ 3 X 10~^G. Nonetheless, even in such cases, the age of the supernova could be thousands 
of years and the enhanced magnetic field is usually much smaller than what is obtained 
by naive cquipartition arguments. In fact, it was clearly anticipated by Meszaros, Rees, & 
Papathanssiou (1994) that the cquipartition argument can at best serve as a broad guide 
to determine the actual value of Bsh, and accordingly, they introduced a completely free 
parameter, A < 1, tagged onto the naively obtained value of B^^: 

B^i^r^Ax lO'^nl^^sX^/'^ G (5) 

Further, it could be possible to apply the basic cquipartition idea at i = to the initial 
FB, and then evaluate the value of the instantaneous Sps or Bsh by using the flux- freezing 
condition. And again, in this case, if we symbolize our ignorance through the free parameter, 
^ < 1, it follows that (Meszaros, Rees & Papathanassiou 1994) 

Bsh - m''^E,,"\Tvl G (6) 
Thus, in this paper, we would attempt to invoke a known feature of relativistic coUisionless 
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MHD shocks to obtain physically significant upper limits on Bgh- 



2. Shock Dynamics 

Following the work on relativistic strong shock jump conditions by Taub (1949) and 
Blandford & McKee (1976) it was discussed in M98 that the ratio of the comoving particle 
densities in the region (1) and (2) is 

r4 ^ r27i2 + 1 . . 

where r2 ~ 4/3 is the effective ratio of the specific heat of (2) and 712 is the bulk LF of (2) 
with respect to (1). Therefore, we have 



n 



2 ~ (4712 + 3)ni 4712711 p=i 2v/27Fni (8) 



where we have used the fact that the maximum value of the LF of the (forward) shocked 
fiuid with respect to the FB is only V2 (M98) : 

^ < 712 < 7f (9) 

The Fig.l actually represents a lab- frame vision where the shock front Si is moving. 
However, shock dynamics is often better studied in the rest frame of Si (*) where, the 
region (1) (upstream) would be seen to gushing in with a large bulk LF, 71* = 7*. On 
the other hand, the region (2) (downstream) would be seen to run away with small LF 
72* ^ (M98). The shock jump conditions also show that the proper relativistic internal 
energy density of the shocked fiuid is 

e'2 ~ 7i27T'2C^ ~ mriirf'c? (10) 

which shows that an equipartition argument in region (2) approximately leads to a value of 
Bsh. as is given by eq. (1): 

B'^ = ~ 150773/^ n]''^ G (11) 
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If the plasma is supposed to have infinite conductivity, all comoving electric fields must 
vanish. Then, for a parallel shock where the embedded magnetic field is parallel to the fiow 
direction (x-axis) one obtains the following simple shock jump condition (de Hoffman & 
Teller 1950, henceforth HT50): 

= B[,; = = (12) 

Further, since under Lorentz transformation, the parallel component of magnetic field 
remains unchanged, for a parallel shock, we find 

— -Bix* = -B2X = -B2X* = Bsh (13) 



2.1. Perpendicular shock 

For a perpendicular shock, the magnetic field changes under Lorentz transformation 
and accordingly shock jump conditions are different from their parallel counterpart. If the 
magnetic field lies in the y-direction (Sx — B^ — 0), Lorentz transformation leads to (HT50) 

where 72* ^ \/2 is the LF of the downstream fluid in the shock frame (M98), and all other 
components of magnetic field in any other frame becomes zero. Further for all transverse 
cases, the frozen-in magnetic field line densities become proportional to the respective 
comoving particle densities (HT50): 



B^^r^ 
B[y n[ 

Then eq. (14) and (15) together imply 



(15) 



^ = ^4-4 (16) 

-Dly 7* "1 

Thus, for the transverse case too, the magnetic field jump condition looks similar to the 
corresponding non-relativistic case if quantities are measured in the shock frame. However, 
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actually, there is a strong enhancement of the magnetic field in the downstream and this 
becomes clear if we write the above equation in terms of respective proper fields by using 



Although, the micro-physics of relativistic shocks, in particular, coUisionless MHD 
shocks is poorly understood, atleast in comparison to their non-relativistic counterparts, 
there have been considerable amount of observational, numerical as well as analytical studies 
of relativistic coUisionless shocks. Although, such studies have largely been focussed for 
the perpendicular or transverse shocks, presumably, to harness the shock drift mechanism 
inherent in the v x B term, we will see that some fundamental aspect of such studies may 
be extendable to the the parallel cases too. To be precise, it is found that the following 
parameter 

where is the speed of the upstream fiuid in the shock frame, is very important for 
studying the structure and accelerating properties of the shock (Kennel & Coroniti 1984a, 
1984b, Alsop & Arons 1988, Hoshino et al. 1992). Physically a is the ratio of the upstream 
Poynting energy fiux and the the particle energy density as measured in the shock frame. 
Stated this way, it may appear that a may not have any relevance for the parallel case. 
However, by noting that B^^ = 'y*B[^ = ^^B[ and that the shock frame upstream particle 
density ni* = 7*ni, one can reformulate eq. (18) as 



eq. (8): 



^2 fa 47i2-Bi fa 47*^1 2V27F-Bi ^ V2r)B[ 



(17) 



3. 



Relativistic MHD shocks 




(19) 



-9- 



One can clearly see that the above definition of a is Lorentz covariant and of a general 



Therefore, the foregoing definition of a becomes relevant for all cases. Also, note that, in 
contrast to non-relativistic shocks, for all practical purposes, ultrarelativistic shocks are 
bound to be tranverse ones because, the shock may be considered to be parallel only if 
^Bi < 7^^, where, is the angle between the upstream (proper) magnetic field and the 
flow direction as measured in the upstream proper frame (Begelman & Kirk 1990, Hoshino 
et al. 1992). In the present case, the upstream proper frame is the ISM and there is no 
reason that the ISM field (which is in any case distorted over certain length scale) should 
be practically perfectly aligned with the GRB fiow direction. Therefore, irrespective of the 
Lorentz covariant definition of cr, the idealized discussion of a perpendicular relativistic 
shock becomes quite relevant in the present case. 

Note that, in case one is performing numerical experiments to study the formation of 
shocks out of relativistic fiows, one must consider a wide range of cr including those > 1 
(Langdon, Arons, & Max 1988). But, it should be realized that a value of cr > 1 necessarily 
means that the upstream flow has enough internal energy to support current systems 
whose energy density exceeds the rest mass energy density. In other words, such numerical 
experiments correspond to an upstream which is relativistically, "hot". For most of the 
realistic astrophysical situations, theoretical and observational arguments suggest that the 
value of a ^ 1 (Piddington 1957, Rees & Gunn 1974, Kennel & Coroniti 1984a,b., Hoshino 
et al. 1992). And, in any case, for the cold upstream region, which is the appropriate for 
the present case and most astrophysical situations, an absolute upper limit is cr < 1. If we, 
somewhat naively, use this upper limit on cr for a parallel shock, we obtain 



nature : 



upstream proper magnetic energy density 



(20) 



cr = 



upstream rest mass energy density 




(21) 
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And this leads to an absolute upper limit for a supposed parallel relativistic shock is 

^B'^^B[^ 0.2ny^ G (22) 

The basic reason that for the parallel case, we are considering the above mentioned upper 
limit on a is that, we are not aware of any theoretical or numerical study of relativistic 
parallel shocks which may tighten this constraint further. Coming back to the fairly well 
studied case of perpendicular relativistic shocks, we may adopt a generous upper limit on 
the value of cr ~ 0.002 by simply adopting the value corresponding to that of Crab (Rees & 
Gunn 1974, Kennel & Coroniti 1984a,b): 

^icrab ^ 9 X 10-^nf G (23) 

Or, 

= ~ 4712^;,,,^ ~ 15 77371^' G (24) 

Given this generous upper-limit on a, depending on the actual obliquity of a given shock, 
the actual value of Bsh will vary between what is shown by eq. (22) and (24) with more 
hkelyhood of assuming the latter value. This condition that we must have cr <S 1 (atleast 
for perpendicular shocks) in order to have a strong and accelerating shock is somewhat 
akin to the condition that for non-relativistic diffusive shock acceleration the Alfven Mach 
number Ma ^ 1 (Begelman & Kirk 1990). This condition means that the Alfven speed 
must be << than the flow speed to ensure that MHD scattering centers scatter the test 
particles vigourously and isotropically. This condition also implies that the magnetic energy 
density is negligible compared to the kinetic energy flux: 

^ << Irmn'ovl (25) 

For nonrclativistic perpendicular shocks too, similar conditions are necessary to ensure that 
the shock is "strong" (Kundt & Krotschek 1982, Leroy et al. 1982, Appl & Camenzind 
1988). 



4. Discussion 



Having obtained this generous upper-limit let us now ponder how justified we are in 
adopting a value of u appropriate for Crab. Remember that the Crab shock is practically 
a standing shock and the upstream region is not the bare ISM. On the other hand, the 
upstream comprises the plasma ejected by the Crab pulsar during its life time of nearly 
thousand years. It is very much likely that, in the past, the value of a for Crab was much 
lower, and the present value has been slowly built up over these thousand years. In contrast 
the case of the GRB blast wave ploughing through the bare ISM is quite difi^erent, and, it is 
highly improbable that, within a lab frame time scale of few seconds or less, the upstream 
medium ahead of the shock front (as seen by in the lab frame) can raise its magnetic 
field from a value of ~ 3 x 10"^G to ~ 9 x 10~^G. We feel that, instead, the following 
scenario is more plausible: The value of B[ probably remains close to its unperturbed 
value ~ 3 X 10~^G; however the shock could be near perpendicular resulting in a value of 
-^sh ~ ^luB'i ~ \/2r]B[ ~ 410^'^?73G. It is also probable, there may not be any stable shock 
formation at all invalidating the rigid relations between B[ and B2 employed so far, and 
on the other hand there may be instantaneous spikes in the downstream magnetic field 
(Langdon, Arons, & Max 1988) : 

%^^[l + {2/a)f' (26) 

and which may erratically raise the shocked field to a Smax ~ 0.1 — IG. However, formation 
of a shock-like discontinuity, either steady or fluctuating requires that the ambient medium 
should be such that the leading particles of the FB, i.e., the piston driving the shock, either 
individually or collectively impart significant amount of their momentum on the ambient 
medium. We endeavoured to examine this critical but usually overlooked problem in M96. 
Noting the similarity between the present problem and the one involving propagation 
of high energy cosmic rays in the ISM, and also recalling that phenomenologically and 
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observationally obtained parameter, the spatial diffusion coefficient describes the entire 
collective interaction of the cosmic rays and the ISM, we found, that, it is implausible that 
the FB can produce any thing akin to a GRB blast wave in the ISM (M96). This is so for 
the simple reason that the BF-ISM interaction time scale estimated in this way could be as 
large as ~ lO^s ! 

Probably, we may also examine here whether, the electrons or the left over pairs of the 
FB, rather than the protons, might carry out the job of imparting energy onto the ambient 
medium (B. Paczynski, private communication). For a given saturation bulk LF of the FB 
(7f ~ 77/2), the electrons or positrons will be less energetic by a factor of ~ 2000. And, if 
we are assuming the Bohm hmit of the diffusion coefficient, the deflection length as well as 
the lab frame deflection time will accordingly be smaller by a factor of ~ 2000. If the numer 
of leptons in the FB could overwhelm the protons by a factor of ~ 2000, there would have 
been equal amount of energy residing in protons and leptons, and, the foregoing reduction 
in time scale would have really meant that the FB-ISM interaction time would be reduced 
by the same factor. Unfortunately, this is not the case, the number of leptons per proton at 
the saturation stage of the optimally baryon polluted FB is indeed ^ 1. This means that 
the eventual FB-ISM interaction time scale should be what was obtained in M96, i.e., 

~ 3 X 10^3 s (27) 

unless it is found that the lepton-ion energy transfer time scale is smaller than this above 
time scale. It may be probable that at best there may be soliton like discontinuites. At any 
rate, we do not think, we have the final answer of such questions at this moment, and we 
realize that, we are simply attempting to understand various aspects subject to our (present 
author's) available knowledge. 

The upshot of this discussion is that the original GRB is likely to be produced either 
by internal collisions within the FB as has long been suspected (Paczynski & Xu, Rees 
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& Meszaros 1994, Papathanssiou & Meszaros 1996) or if the original FB is propagating 
within a medium which has a modest baryonic mass but is dense enough to absorb the FB 
momentum either by binary coUisions or by coUisionless MHD process. Finally, reverting 
back to the original mandate of this paper, in the case of a supposed GRB blast wave 
propagating in the ISM and even ignoring the disturbing possibility that no shock like 
discontinuity may be formed on the GRB time scale, we feel that, the maximum value of 
the magnetic field in the shocked fluid may not exceed ~ IG. We would like to emphasize 
here the fact that this conclusion does not at all imply that the FB can not excite a blast 
wave (initially relativistic and then non-relativistic) on a much larger time scale of days 
or weeks and which is necessary for explaining the GRB after glow in various low energy 
brackets. 
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Figure Caption 

Figure 1: Sketch of the FB-shock configurations, for details, see text. 



This manuscript was prepared with the AAS macros v4.0. 



4 

> 

FB 



2 

> 



1 

ISM 

S1 



